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ABSTRACT 

A detailed light curve analysis is presented for a first-overtone RR Lyrae star, MACHO* 
J050918.712-695015.31, based on MACHO and OGLE-III observations. As a foreground 
object of the Large Magellanic Cloud, it gives an extraordinary opportunity to study an almost 
continuous, 12-year long dataset of a relatively bright (V « 15™0) RRc star with rapid period 
change. Cyclic period modulation is suggested by the O — C method, where the cycle length 
is about 8 years. With the available unique dataset, we could draw strong limits on other 
light curve changes that may be associated to the period modulation. We could exclude both 
multiple periodicity and amplitude modulation unambiguously. Any theoretical model should 
reproduce the observed lack of photometric modulations. Simple arguments are also given for 
possible hydromagnetic effects. 
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1 INTRODUCTION 

Period changes in RR Lyrae stars have been an intriguing issue 
since the discovery of this phenomenon (Bailey 1913, Prager 1939). 
Numerous studies have been carried out for RR Lyrae stars in glob- 
ular clusters (e.g., Szeidl 1975, Wehlau et al. 1992, Rathbun & 
Smith 1997, Jurcsik et al. 2001) and in the galactic field (see a 
collection of field stars in Firmanyuk 1976, 1982, and a recent 
re-analysis of a particular object by Jurcsik et al. 2002). Stellar 
evolution is expected to cause very slow continuous changes in 
period (Lee 1991, Cox 1998), but the observations tend to reveal 
abrupt and/or cyclic changes that are too fast to be explained in this 
way. The Blazhko effect (Blazhko 1907), which refers to ampli- 
tude modulation in RR Lyrae stars, does not yet have a generally 
accepted explanation (Chadid et al. 2004) and is also known to be 
associated with period modulations. To make the overall picture 
even more complicated, Paparo et al. (1998) found opposite period 
changes in the modes of two double-mode RR Lyrae stars in M15. 
A number of theoretical studies have attempted to explain period 
changes in RR Lyrae stars with various mechanisms, such as mix- 
ing events in the semiconvective zone (Sweigart & Renzini 1979). 
Good reviews can be found in Rathbun and Smith (1997) and Smith 
(1995, 1997). Understanding these period variations remain an im- 
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portant unsolved problem in the stellar pulsation theory (Alcock et 
al. 2000). 

Globular clusters that are rich in RR Lyrae variables offer a 
good opportunity to study period changes. For example, Rathbun & 
Smith (1997) studied the period changes of RR Lyrae variables in 
seven globular clusters by collecting all data in the literature over 50 
years. They found that the rates of period changes of fundamental- 
mode RR Lyrae (RRab) stars are usually larger than those of first- 
overtone (RRc) and double-mode (RRd) variables. However, they 
concluded that current theories did not explain the irregular period 
changes. More recently, Jurcsik et al. (2001) presented an analysis 
of the RR Lyrae variables of uo Centauri. They studied O — C di- 
agrams of altogether 126 RRab, RRc and RR Lyrae-like stars and 
found that the period changes of most RRab were in good agree- 
ment with the evolutionary model predictions, but those of the first- 
overtone RRc-type stars showed a much more complex, irregular 
behaviour (Jurcsik et al. 2001). 

Despite the great value of cluster RR Lyrae stars, no one has 
ever followed a cluster continuously with the same instrument for 
many years or decades. The essentially random samples over many 
decades with various instruments do not provide a homogeneous 
view of the period changes of RR Lyrae variables. The situation 
changed in the early 1990s, when microlensing projects started 
monitoring the sky towards the Large and Small Magellanic Cloud 
and the Galactic Bulge, providing continuous photometric obser- 
vations of millions of stars during several years. Of these, MA- 
CHO and OGLE projects had the best coverage. These data present 
a great opportunity for investigating period changes of RR Lyrae 
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stars. Alcock et al. (2000, 2004) analysed MACHO light curves for 
more than 1000 RRc stars and found that 10% showed strong pe- 
riod changes, whose strength and patterns excluded simple evolu- 
tionary explanations. Meanwhile, Soszyhski et al. (2003) compiled 
a more complete catalog of ~ 7600 RR Lyrae stars in the LMC 
based on the OGLE-II observations. It remains for a future study 
to combine these two sources of data for examining period changes 
over a longer time base. 

Although a large amount of data is available for the RR Lyrae 
stars in the Large Magellanic Cloud, there is a major drawback 
with these stars: they are very faint at the distance of the LMC 
(V « 19^0). Hence, the signal-to-noise ratio (S/N) of the MACHO 
observations for the LMC RR Lyrae stars is very low (Alcock et 
al. 2000) and neither light curve parameters nor their changes can 
be determined accurately. Although OGLE-II RR Lyrae observa- 
tions have higher S/N ratio due to better seeing and more advanced 
photometric reductions technique, they were carried out over a 
shorter time span and with only a single filter. These problems 
could be circumvented if we found suitably bright foreground stars 
in the MACHO and OGLE databases. During an analysis of more 
than 6000 MACHO variable stars (Derekas et al. in prep.), we no- 
ticed MACHO* J050918.712-695015.31 (hereafter M0509-69), 
an RR Lyrae star with intriguing period changes. At magnitude 
V « 15™0, period P = 0.328d and amplitude A hluc « (T5, 
this star offers a unique opportunity, because of its high S/N (~ 
80 - 120 in the MACHO data and ~ 120 - 130 in the OGLE- 
III data) light curve, for analysing the correlation between period 
changes and light curve shape variations. This is expected to tell us 
something about the physical mechanism causing the modulations. 
We note that M0509— 69 was also included in the study of Alcock 
et al. (1997), who discussed 20 foreground RR Lyrae stars towards 
the LMC and noticed the changing period of this star. 

The main aim of this work is to demonstrate the period mod- 
ulation of M0509— 69, and show that it is not accompanied by any 
variations in amplitude or light curve shape. Our study exploits the 
homogeneity, the continuity and the 12-year long time span of the 
combined MACHO and OGLE-III CCD observations. In Sect. 2 we 
describe the data and their basic analysis. We show several pieces 
of evidence for pure period modulation of M0509— 69 in Sect. 3. 
Finally, we briefly discuss possible mechanisms in Sect. 4. 



among the first-overtone RR Lyraes (Poretti 2001), the light curve 
showed enormous phase modulations. 

In light of the strange behavior of M0509-69 in the MACHO 
data, we decided to supplement them with OGLE observations. Un- 
fortunately, this object is located outside the fields covered during 
the second phase of the OGLE project (Udalski et al. 1997) so that 
OGLE-II photometry from 1997-2000 is not available. However, 
this field was covered in the OGLE-III phase, which started in June 
2001, and the star is continuously monitored up to now. Thus the 
gap between the MACHO and OGLE-III observations is only 621 
days. The OGLE observations were obtained with the 1.3 m War- 
saw telescope located at Las Campanas Observatory, Chile, which 
is operated by the Carnegie Institution of Washington, equipped 
with 8192 x 8192 pixel mosaic camera. Photometry was derived 
with the standard OGLE data pipeline based on DIA image sub- 
traction technique (Udalski 2003). Observations were done through 
/-band filter. 

As the first step of the analysis of M0509— 69, we searched 
for the period using the Phase Dispersion Minimization method 
(Stellingwerf 1978), then we refined that period using the String 
Length method (Lafler & Kinman 1965, Clarke 2002). The phase 
diagram with the finally adopted period (Sect. 3.2) is plotted in 
Fig.Q Both the phase diagram and random checks of short sub- 
sets suggested a strong period modulation, so we therefore deter- 
mined epochs of maximum light to construct the classical O — C 
diagram. To this, we divided the whole dataset into 50-day long 
subsets, which did not show any noticeable phase shift. Then we 
phased each subset with the period given by the String Length 
method (0.32806 days). When there were too few points within 50 
days, then we used 100-day long subsets. The best of these phase 
diagrams was fitted with a fourth-order Fourier-polynomial. This 
polynomial was used as a master curve - allowing both vertical 
and horizontal shifts to fit every other phase diagram. The same 
procedure was applied for the OGLE-III data. From the resulting 
phase shifts, the period and the epochs of each phase diagram, we 
calculated the times of maximum light, listed in Table Q Typical 
uncertainties are about 0.001-0.0015 days. 

Finally, we note that downloadable MACHO data are in Mod- 
ified Julian Date, so that we had to add 0.5 days to all epochs of 
maximum. 



2 INITIAL DATA ANALYSIS 

Our attention was drawn to this star during an on-going study of 
MACHO variable stars, that are classified as eclipsing binaries in 
the on-line MACHO Variable Star Catalog. The MACHO observa- 
tions were obtained between June 1992 and January 2000, with the 
50 inch telescope at Mt. Stromlo Observatory using a two-channel 
system (Hart et al. 1996). Two-colour data in the specifically de- 
signed MACHO blue (B M ) and MACHO red (i? M ) bands are pub- 
licly available at the MACHO website . Stars have been classified 
by an automatic software and we downloaded individually every 
target, 6837 in total. However, it became obvious very quickly that 
the classification was not perfect, as large fraction of "eclipsing bi- 
naries" was found to be Cepheids, RR Lyrae stars or long-period 
variables. We therefore decided to re-classify all stars and during 
this procedure we noted the peculiar light curve of M0509— 69. 
Whereas the period and amplitude unambiguously place this star 

1 |http : / /wwwmacho .monaster . ca/| 



3 RESULTS 

3.1 Multiple periodicity 

We examined the possibility of multiple periodicity in the MA- 
CHO blue-band light curve of M0509-69. We chose this set be- 
cause of the higher amplitude, thus better S/N. The frequency anal- 
ysis was carried out by standard Fourier-analysis using Period98 of 
Sperl (1998). The calculated amplitude spectrum is shown in the 
top panel of Fig. [2] The main peak at /o = 3.048222d _1 is in 
very good agreement with the result by the String Length method 
(3. 048223d -1 ). Usual prewhitening steps resulted in the same be- 
haviour as the one found for RR1 -PC-type stars by Alcock et al. 
(2000): the prewhitened frequency spectrum contained significant 
remnant power very close to the main component. As was shown 
in Alcock et al. (2000), this can be attributed to long-term pe- 
riod and/or amplitude change. Since the procedure of determining 
epochs of maximum proved the stability of the amplitude (i.e. ev- 
ery phase diagram could be fitted very well with the master curve), 
we checked for the presence of multiple periodicity in two different 
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Figure 1. The phase diagram for M0509-69 in three colours: MACHO blue, MACHO red and OGLE. 



Table 1. Times of maximum for M0509-69. 
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ways. Firstly, we did regular prewhitening in the best three 50-day 
long subsets to check the residual frequency spectra for any signif- 
icant peak unrelated to the pulsational frequency. For that reason, 
we fitted and subtracted not only /o, but 2/o, 3/o and 4/o from the 
subsets. None of the residual spectra contained a significant peak 
other than 5/o, after which we could not find anything else. 
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Figure 2. The Fourier spectra of the complete blue colour dataset (top 
panel) and the residuals, after removal of the main period, harmonics and 
aliases (bottom panel). The insert shows the window function on same fre- 
quency scale. 



Secondly, we subtracted individually 
polynomials from every subset in form of 



A + A * 



sin (2irift + < 



fitted Fourier- 



(1) 



where the amplitudes and phases were fitted with Period 98. After 
subtracting each polynomial, we re-united the residual light curves 
and calculated the frequency spectrum of the whole set. As ex- 
pected from the asymmetric shape of the mean light curve, we de- 
tected only 5/ and 6/ with amplitudes of 0^0064 and 0^0045. 
After subtracting these two frequencies, the S/N ratios (Breger et 
al. 1993) of the next three peaks were about 4, close to the limit of 
significance. We could identify these frequencies as close remnants 
of the primary peak and its harmonics or daily aliases. Finally, we 
arrived at pure white noise (bottom panel in Fig.|2j, confirming the 
monoperiodic nature of M0509— 69. Similar results were given by 
the red data too, but their S/N is lower than that of the blue data, 
which we preferred in our analysis. 
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3.2 The O C diagram 

Having excluded multiple periodicity, the O — C diagram was made 
from the times of maxima in Table ^Initially, we used the period 
given by the String Length method, with the following ephemeris: 
-ffJAnax = 2449197.4241 + 0.32806 x E. It showed very clear 
cyclic period changes without colour dependence (i.e. the blue and 
the red O — C points drew the same curve). 

There is a well-known problem with the O — C method: al- 
though it is a very effective tool to detect tiny period changes, the 
shape of the diagram strongly depends on the period used in the cal- 
culation. Moreover, it can be discontinuous after long periods, so 
that we analysed the O — C diagram using the continuous method 
suggested by Kalimeris et al. (1994). In this approach, both the pe- 
riod and its rate of change are continuous functions of time. Rather 
than using various assumptions on period change (e.g. assuming 
constant change of period and thus fitting parabola to the O — C\ 
or fitting straight lines to selected parts of the O — C), we deter- 
mined the "instantaneous" periods from every neighbouring pair of 
subsets, simply as the local derivative of the O — C plot plus the 
period used in the ephemeris (P e ). The difference between the "in- 
stantaneous" period and P e is plotted in the bottom panel of Fig.[3] 
As can be seen here, the period change is enormous: the difference 
between the longest and the shortest period is about 12 seconds; in 
other words, the relative period change is ±2 ■ 10~ 4 ! 

The use of the "instantaneous" period has a further advan- 
tage: this way we can combine MACHO and OGLE-III data with- 
out involving systematic phase shift between the Bm and I bands. 
We checked several foreground RR Lyraes with simultaneous MA- 
CHO and OGLE-II observations (Zebrun et al. 2001) and found 
that I-band maxima are usually within 0.01—0.03 in phase relative 
to i?M-band maxima. With the "instantaneous" period, even this 
small effect is eliminated. 

Finally, we corrected the period by assuming that the time- 
span of the combined MACHO+OGLE-III data is long enough to 
average out the period modulation. In other words, we searched that 
period which resulted in zero mean period difference in the bottom 
panel of Fig. [5] The resulting ephemeris is the following: 



HJD„ 



2449196.9241 + 0.328038(2) x E, 



(2) 



and the finally adopted O — C diagram in the top panel of Fig. [5] 
was also calculated with this ephemeris. 



3.3 Amplitude modulation 

Some RR Lyrae stars show cycle-to-cycle changes in the shape 
of their light curves, which can be caused either by double-mode 
pulsation (in RRd stars) or the still-puzzling Blazhko-effect. Since 
double-mode pulsation was ruled out in Sect. 3.1, what remains to 
be checked is the possibility of amplitude modulation. 

Although the master curve fitted very well all subsets, 
we wanted to examine quantitatively the possibility of light 
curve shape change. Therefore, we fitted fourth-order Fourier- 
polynomials to each subset using equation (1) at the primary fre- 
quency and its harmonics. Then we measured the light curve shapes 
with the Fourier parameters (Simon and Teays 1982), of which we 
show particular results for R21 = A2/A1 and R31 = A3/A1 
(these are the most accurate parameters). Fig. [4] shows the relative 
variations of R21 and R31. The error bars are quite large (because 
the light curve is not too asymmetric, so that higher harmonics are 
fairly weak), but it is evident that there were no measurable shape 
changes. 
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Figure 3. Top panel: the O — C diagram of M0509— 69 from Bm and / 
data. Bottom panel: the "instantaneous" period. The error bars are the 3<r 
errors of the fits. 
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Figure 4. The relative variations of R21 and R31 . 

Since individual subsets rarely contained enough data points 
for determining accurate Fourier-parameters, we also examined the 
correlation between the period change and the light curve shape 
using five subsets around the minimum and five around maximum 
of the "instantaneous" period. Then we corrected these phase dia- 
grams for any phase shifts to form high-quality, continuous phase 
diagrams representing the two extrema of the period. Then we fit- 
ted sixth-order Fourier polynomials to both phase diagrams and the 
resulting parameters are shown in Table |2| As can be seen, none 
of the Fourier-coefficients differs significantly for the two cases. 
This means that we could not detect any light curve change associ- 
ated with the period variation. This conclusion is also supported by 
the phase diagram of the whole MACHO and OGLE-III datasets 
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Figure 5. The phase diagrams after eliminating the phase shifts between all subsets. 



Table 2. The coefficients of the sixth-order Fourier polynomials around the 
minimum and the maximum of the period. 
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-0.313 ±0.476 
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(Fig. H}, where we shifted every subset to match the master curve 
in phase. The well-defined shapes show no evidence for any addi- 
tional variation in light curve shape. 



4 DISCUSSION 

It is obvious from Fig.|5|that the period change seems to be cyclic, 
which implies there must be some (quasi)periodic mechanism af- 
fecting the pulsation. The shape of the curve in the bottom panel of 
Fig.[5]is largely independent of the period used for the O — C dia- 
gram (provided that it is close to the real mean period), so that the 
observed change cannot be explained by some numerical artifact 
caused by a wrong period in the ephemeris. 

As a first possibility, we consider light-time effect (LITE) due 
to orbital motion in a binary system. Although it was suggested 
for several RR Lyrae stars that their cyclic O — C diagrams were 
caused by LITE, so far this was supported by spectroscopic ob- 
servations only for TU UMa (Wade et al. 1999). For M0509-69, 
the half-amplitude of the O - C in Fig. is A(O-C)=0.09 d, 
while the time separation of the two minimum is 2900 d. Assuming 
that A(0 — C) = as >™\ the estimated semi-major axis is about 
15.5/ sin i AU. Given P or b=2900 d « 7.95 yr; the minimum mass 
of the system would be almost 60Mq . 

Although there are models that predict the existence of de- 
tached black hole binaries with Mbh = 10 — 20 Mq and M sec = 
0.3 — 1 Mq (see, e.g., Table 1 in Podsiadlowski et al. 2003), the sec- 



ondary components in these theoretical systems are not expected to 
follow the evolution of "normal" stars with the same masses. This 
is because they are descendants of more massive progenitor stars 
that transferred most of their masses to the black hole primaries 
(but models also predict that they can evolve through the classi- 
cal instability strip, even becoming "RR Lyrae-like" variables). In 
principle, we cannot exclude this possibility until it is checked by 
spectroscopy (any sign of orbital motion should be relatively eas- 
ily detectable spectroscopically even at V = 15 mag), although 
6OM0 seems extraordinary large. Also, the high occurrence rate of 
strong period changes in RR Lyrae stars argues against this exotic 
explanation. 

Stellar evolution can be safely excluded, since both the 
strength and the cyclic nature of the period change exceed theo- 
retical expectations by orders of magnitudes (Smith 1995). 

Another interesting possibility is related to hypothetic hydro- 
magnetism. Stothers (1980) argued that RR Lyrae stars may be 
considered as fair analogues of the Sun. He proposed that radius 
changes may be driven by some sort of magnetic activity, which in 
turn would cause observable period changes. The correlation be- 
tween the frequency changes and the magnetic variations was con- 
firmed by Howe et al. (2002), who performed a detailed study on 
the solar cycle frequency shifts in global p-modes of the Sun. They 
found that the latitudinal distributions of the frequency shifts shows 
close temporal and spatial correlation with the unsigned surface 
magnetic flux. 

In our case we can calculate the radius change of M0509— 69 
from the period-density relation, i.e. P^fp = Q. Substituting the 
values determined from the bottom panel in Fig. [3] it follows: 



AR 
R 



AP 



28342 



1.4 • 10" 



(3) 



Note that this equation is an approximation that omits a term that 
contains the magnetic energy (see Stothers 1980). 

Very recently, Noel (2004) presented an analysis of 13 years 
of apparent radius measurements of the Sun. Interestingly, he found 
a correlation between the solar radius variation and the sunspot 
numbers, related to the solar magnetic activity. Based on the data 
in his fig. 5, the relative change of the solar radius was about 
AR/R ~ 4 ■ 10~ 4 . Assuming that RR Lyrae stars can show 
solar-like magnetic activity cycles, these two similar relative ra- 
dius changes seem to support the hypothetic hydromagnetic period 
change. 

We know little about the presence of magnetic fields in 
RR Lyrae stars. Most recently, Chadid et al. (2004) presented high- 
precision longitudinal magnetic field measurements of RR Lyrae 
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itself taken over 4 years. Their results provided no evidence for a 
strong magnetic field in the photosphere of this star, which was of- 
ten postulated in various explanations for the Blazhko effect (see, 
e.g., Jurcsik et al. 2002). If the period change in M0509-69 is in- 
deed caused by solar-type magnetic activity, this field might be 
quite difficult to detect observationally. In fact, the only Blazhko- 
related phenomena that have been detected so far seem to be small 
line-profile variations (possibly related to non-linear pulsations; 
Chadid et al. 1999, Chadid 2000) and even these are only detectable 
with high S/N high-resolution spectra. 

Although the effects of hydromagnetism on pulsating stars are 
not well understood, we conclude that the Stothers-mechanism can 
perhaps explain the pure period modulation of M0509— 69. How- 
ever, new theoretical investigations are necessary to interpret this 
observed behaviour. Detailed models should also address the ab- 
sence of any light curve shape change and multiple periodicity. 
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